Perturbative description of nuclear double beta decay transitions 
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A consistent treatment of intrinsic and collective coordi- 
nates is applied to the calculation of matrix elements describ- 
ing nuclear double beta decay transitions. The method, which 
was developed for the case of nuclear rotations, is adapted to 
include isospin and number of particles degrees of freedom. 
It is shown that the uncertainties found in most models, in 
dealing with these decay modes, are largely due to the mix- 
ing of physical and spurious effects in the treatment of isospin 
dependent interactions. 
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One can hardly overestimate the importance of the 
double beta decay as a process explicitly linking the 
physics of neutrinos with the nuclear structure Jl]-||] . Un- 
fortunately, the reliability of the theoretical predictions 
has been hampered by unstabilities in the many-body 
BCS + RPA - type of treatments that have been applied 
during the last decade ftfl- An alternative approach 
based on group theoretical methods has confirmed the 
existence of a zero-energy state for certain values of the 
strength of the proton-neutron, particle-particle, effec- 
tive interaction j|J^] . The appearance of such a state has 
been interpreted as a signature of a phase transition [pi . 

Here we take an alternative point of view based on 
the fact that the zero-energy state is a consequence of 
the breakdown of the isospin symmetry implicit in the 
(separate) neutron and proton BCS solutions 

As similar to the case of deformed nuclei, the sym- 
metry may be restored in the laboratory frame through 
the introduction of collective coordinates. However, as 
different from previous cases dealing with collective co- 
ordinates, we must use here an interaction which does not 
conserve isospin, in order to obtain allowed matrix ele- 
ments between T — > (T — 2) states: a central many-body 
problem that must be solved in double beta decay calcu- 
lations is to disentangle unphysical isospin violations in- 
troduced by the theoretical treatment from physical effects 
due to the interaction. 

In this letter we present the formalism for the case of 
particles moving in a single j-shell and coupled through a 
monopole pairing force. The hamiltonian is not necessar- 
ily isoscalar. This simple case involves all the complica- 
tions associated with the collective treatment. Moreover, 
the predictions may be compared with those of exact cal- 
culations p] of nuclear double beta decay transitions of 



the Fermi- type. 

We define the operators 



pm nm ' nm pm / 



m>0 m>0 

T A = \{t p + T n ) ; r = -(r p - t„) ; [tj, t{\ = tq, 



To — - (T p — T n ) ; Pi, Ti] — — T , 
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T A = \{T P + T n ) 

eA — e p + e„ ; e = e p - e„. (1) 

where t v and e v are the number operators and single- 
particle energies, respectively, and v = p,n, 1 = — 1. 
The T v ,T±i are the generators of collective rotations in 
gauge-and isospace. 

The collective treatment appropriate for an isospin 
conserving pairing interaction was introduced in refs. 
Jl2| , 13 1, The basic set of states associated with 

the collective sector may be labeled by the four quantum 
numbers \Ta] T, to, k >, where Ta is the total number of 
pairs of particles. We substitute M, T (the isospin pro- 
jections in the laboratory and intrinsic frames) by the 
quantum numbers m = \{T + M) and k = |(T + To), 
respectively. We focus on states such that m « T and 
k = 0. Hereof we drop the labels Ta, k from the collective 
states. 

The introduction of collective degrees of freedom is 
compensated through the appearance of the constraints 



t z -T z = 0; (z = n,p,±l), 



(2) 



which express the fact that we can rotate the intrinsic sys- 
tem in one direction or the body in the opposite one with- 
out altering the physical situation jn|. Physical states 
should be annihilated by the four constraints and physi- 
cal operators should commute with them. 

Unphysical violations of the isospin symmetry take 
place in the intrinsic frame, which may be defined, for in- 
stance, by the condition S± = 0, where the bar denotes 
the g.s. expectation value. This condition is precisely 
satisfied by performing a separate Bogoliubov transfor- 
mation for protons and neutrons. 

Physical isotensor operators have to be transformed 
from the laboratory frame to the intrinsic frame. In the 
case of the single-particle and pairing hamiltonians we 
have 



= e A T A + eoiD^ro + D^n + D^n) 
~ eATA + e (r + d + d + t^titj) 



(3) 
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S+S±)] 
(6) 



-g 2 [(S+S p + S+Sn) + - S n Sp([3- 4 d 2 + (3 A d+ 2 ) 



Jp Up 
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<S 2 + S 2 n )(^d+d + -L n 7i) - ^S n Sp(r 2 + T 2 )], 

plus null terms which are proportional to the constraints 
(|J). We have kept only the lowest order terms in an 
expansion in powers of T _1 , assuming 0(S V ) — T and 
0(g v ) — T _1 . Use has been made of the relation S± 

~ — ^{S n T\ + SpTf). The have been expressed 

by means of Marshalek's generalization of the Holstein- 
Primakoff representation |l6|] , such that the operators 
f3,d + ,d satisfy the relations [d, d + ] — 1 and 



(P) 2t \T,m > =\T + t,m>, 



(d + Y\T,m> = 



(m + t)! 



T,m + t> 



(7) 



It is easy to verify that the four components of the hamil- 
tonian (||)-(^) commute with the constraints (||) and are 
therefore physical operators. 

The contributions in (^)-@ that are independent 
of the operators f3, d + , d, t±\ add together to the two 
(proton-proton and neutron-neutron) pairing hamiltoni- 
ans in a single j-shell, namely 



9v S^ S v 



(8) 



with interaction strengths g p = 90 + 91 + 92 and g n — 
9o — gi + 52 • These hamiltonians are separately treated 
within the BCS approximation. Lagrange multiplier 
terms —X v (t v — T v ) are added before solving the BCS 
equations. This treatment yields the independent quasi- 
particle energy terms E v v v = hQ,g v u v , where v v is the 
quasi-particle number operator for i>-nucleons and f2 is 
half the value of the shell degeneracy. Within the RPA 
we may write 



EpVp + E n v n 



rp2 [SpS n {T~i + r x 



-(T 2 



s 2 )Tin 



(9) 



Adding these contributions to the remaining terms 
in (0)-@ (i.e., to the ones depending on the operators 
/?, d + , d, r±i ), one obtains for the proton- neutron sector 
of the spectrum (to leading order in T -1 ) 



H, 



(lab) _ 



(sp+p) 



= H + u d d + d + H 2 , 



(10) 



H -— g v t v g v S V i 

Ud = e + j;(S 2 - SI) + ^-(S 2 + §1), 



< T- 2,m- 2\H 2 \T : m >-- 



3,9s 
T 



S p S n yJm{m - 1), 



plus null terms. 

To leading order, the isospin operators in (||) - (H) and 
in (^) have a boson structure since [tj, (— n)] ~ T and 
tj annihilates the state with fo = —T. This is pre- 
cisely the phonon that yields a zero frequency root for 
isoscalar hamiltonians within a naive RPA @. In the 
present treatment this phonon has disappeared from the 
final physical hamiltonian (10) to become part of the con- 
straints (2). Nevertheless this degree of freedom must be 
taken into account in higher orders of the expansion in 
powers of T _1 , for instance through the BRST procedure 
fL7f , as applied to many-body problems in and to 
the particular case of high angular momentum in |l8| . 

The spectrum of the system is ordered into collective 
bands associated with particular values of the total num- 
ber of particles and the isospin (T < Ta)- The energy of 
the band head is given by the BCS expectation value H. 
The different members of each band are labeled by the 
quantum number m and are separated by the distance 
u>d, which includes the difference between the proton and 
the neutron single-particle energy eo- There is an inter- 
band interaction which mixes different values of T but it 
conserves the projection M in the laboratory frame (cf. 
eq. ©). 

From the point of view of the expansion in powers of 
T _1 , the interband interaction is of the same order (0(1)) 
as the distance between the states that are mixed by 
it. Nevertheless, in the following we continue applying 
perturbation theory and we also request that l^l < 9v 

In the calculations that we report in this paper we as- 
sume 9p = 9n= 9- The excitation energy u>d is displayed 
for the cases j = |, Ta = 5, T = 3, eo = 0.8 MeV, g — 
0.4 MeV and j = i|, T A = 10, T = 4, e = 0.63 MeV, 
g = 0.2 MeV as functions of the ratio (72/5 (upper boxes 
of fig. 1). We predict the exact results for g 2 = and 
very satisfactory ones for the other values, in spite of the 
fact that for these results we have neglected the interband 
interaction. 

The strong current that appears in the weak hamilto- 
nian is proportional to the isospin operator 

/3_ = -\/2rf ab) = -V2(Dl in + D\ t q + Dj T ) 
w V2Td + + null operator. (11) 
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The matrix element of double beta decay transitions, 
which for the present case correspond to pure Fermi tran- 
sitions (cf. ||), is proportional to the product of the two 
matrix elements 

Mi = < T, l\p-\T,0 >« V2T 
M 2 = < (T-2),0|/?_|T,1 > 

2 V / T< (T - 2), 0|iJ 2 |r, 2 > 



H(T,2)-H((T-2),0) 



(12) 



These matrix elements are displayed in the lower boxes of 
fig. 1 for the same parameters as in the upper boxes. The 
expression for the interband matrix element in ( |i"fj| ) docs 
not distinguish whether the r.h.s. should be calculated 
for the initial or the final value of T, since it is valid 
for T » 1. Therefore, the effective interband matrix 
element has been chosen as the geometric average of the 
values obtained for each of the two connected bands. 

Fig. 2 displays Fermi double beta decay matrix ele- 
ments, corresponding to transitions from the initial to 
the final ground states. It has been calculated using the 
expression 



M 2v = 



M X M 2 



(13) 



where the energy released A has been taken to be 0.5 
MeV, as in || . In addition to the exact and perturbative 
values of these matrix elements, we have included in this 
figure the results obtained by using some other approxi- 
mations. The exact result shows the suppression of the 
matrix element around the point where the strength of 
the proton-neutron symmetry breaking interaction ap- 
proaches the value of the fully symmetric interaction. 
This result is reproduced both in the naive QRPA and 
in the perturbative approach. The other approximation 
badly misses this cancellation. A detailed comparison 
between the results of exact, naive QRPA and renormal- 
ized QRPA (RQRPA) calculations can be found in ||. 
It is worth to note that in the perturbative approach the 
corresponding sum rule (Ikeda's sum rule) is exactly ob- 
served. This is not the case of other approaches, like 
in the case of the RQRPA of pTsj] , where the sum rule 
is violated. One can easily understand this failure of 
the RQRPA approach, since it badly mixes-up terms in 
defining the components of the equation of motion. Simi- 
lar conclusions about the validity of the RQRPA method 
are reported in [^j. The perturbative approach, as seen in 
figs. 1 and 2, not only reproduces exact results very sat- 
isfactorily but it also gives some insight about the mech- 
anism responsible for the suppression of the matrix el- 
ements. As found in the calculations, the value of the 
matrix element M 2 depends critically on the strength of 
symmetry breaking terms of the interaction between pro- 
tons and neutrons. The terms are proportional to 172, as 
shown before. On the other hand, the values of Mi are 



not very much dependent on this interaction. Finally, 
it should be observed that the point where the excita- 
tion energy vanishes and the point where the symmetry 
is completely restored are different (cf. fig. 1). This 
result, also obtained in the exact diagonalization of the 
full hamiltonian, cannot be reproduced by other means 
as shown in ||. Further details will be presented in a 
longer publication, in which the extension of the formal- 
ism to include any number of non-degenerate j-shells has 
been performed and an isotensor isospin interaction has 
been included pCfl . 

In conclusion, it is found that a correct treatment of 
collective effects induced by isospin dependent residual 
interactions in a superfluid system is feasible: physical 
effects due to the isospin symmetry-breaking terms in 
the hamiltonian are obtained even in the presence of the 
BCS mean field built upon separate proton and neutron 
pairing interactions. The definition of intrinsic and col- 
lective coordinates and their separation guarantees that 
the isospin symmetry is restored and that spurious con- 
tributions to the wave functions are decoupled from phys- 
ical ones. Particularly, the problem of the unstabilities 
found in the standard proton-neutron QRPA are avoided 
by the explicit elimination of the zero frequency mode 
from the physical spectrum but keeping it in the pertur- 
bative expansion. The appearance of this mode cannot 
be avoided by the inclusion of higher order terms in the 
QRPA expansion or by any other ad-hoc renormalization 
procedure, like the renormalized QRPA of JL9| , once the 
BCS procedure is adopted for the separate treatment of 
proton and neutron pairing correlations. 

The results shown in this letter are very encouraging, 
in spite of the fact that we have not used very large values 
of T. We shall discuss the extension of the formalism 
to cases such that T « 0(S V ) and its application to 
Gamow- Teller transitions in a forthcoming publication. 
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FIG. 2. Matrix elements for Fermi double beta decay tran- 
sitions calculated in several different approximations. The 
meaning of approximations denoted as Quasiparticle Random 
Phase Approximation (QRPA) and Renormalized Quasipar- 
ticle Random Phase Approximation is explained in the text. 
Exact and Perturbative results for the two sets of parameters 
utilized in the calculations are indicated by solid and dotted 
lines, respectively. 



FIG. 1. Excitation Energy and Transition Matrix Ele- 
ments. Exact (solid lines) and Perturbative (dotted lines) re- 
sults for the excitation energy (upper boxes) and transition 
matrix elements Mi and Mi (lower boxes) corresponding to 
the two different sets of parameters (j = 9/2 and j = 19/2) 
discussed in the text. 
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